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Abstract Based on the Sloan Digital Sky Survey DR6 (SDSS) and Millennium Simulation 
(MS) we investigate the alignment between galaxies and large-scale structure. For this pur- 
pose we develop two new statistical tools, namely the alignment correlation function and the 
cos(29)-statistic. The former is a two-dimensional extension of the traditional two-point cor- 
relation function and the latter is related to the ellipticity correlation function used for cosmic 
shear measurements. Both are based on the cross correlation between a sample of galaxies 
with orientations and a reference sample which represents the large-scale structure. We ap- 
ply the new statistics to the SDSS galaxy catalog. The alignment correlation function reveals 
an overabundance of reference galaxies along the major axes of red, luminous (L > L») 
galaxies out to projected separations of 60 /i _1 Mpc. The signal increases with central galaxy 
luminosity. No alignment signal is detected for blue galaxies. The cos(2#)-statistic yields 
very similar results. Starting from a MS semi-analytic galaxy catalog we assign an orienta- 
tion to each red, luminous and central galaxy, based on that of the central region of the host 
halo (with size similar to that of the stellar galaxy). As an alternative we use the orienta- 
tion of the host halo itself. We find a mean projected misalignment between a halo and its 
central region of ~ 25°. The misalignment decreases slightly with increasing luminosity of 
the central galaxy. Using the orientations and luminosities of the semi-analytic galaxies we 
repeat our alignment analysis on mock surveys of the MS. Agreement with the SDSS results 
is good if the central orientations are used. Predictions using the halo orientations as proxies 
for central galaxy orientations overestimate the observed alignment by more than a factor of 
2. Finally, the large volume of the MS allows us to generate a two-dimensional map of the 
alignment correlation function which shows the reference galaxy distribution to be flattened 
parallel to the orientations of red luminous galaxies with axis ratios of ~ 0.5 and ^0.75 for 
halo and central orientations, respectively. These ratios are almost independent of scale out 
to60/i _1 Mpc. 

Key words: dark matter halos: clustering - galaxies: large-scale structure of Universe - cos- 
mology: theory - dark matter 



1 INTRODUCTION 

Recent large redshift surveys, lik e the 2dF Galaxy Redshift Survey (2dFGRS. IColless et al.ll200ll) and the 
Sloan Digital Sky Survey rSDSS. lYork etal]l2000b . allow the cosmic large-scale density field to be traced 
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with unprecedented accuracy. Different structures can reliably be classified as groups, filaments, walls and 
voids and match well with the patterns seen in N-body simulations. These structures induce large-scale 
tidal fields which in turn cause large-scale correlations in the orientations of massive dark matter halos (cf. 
Bon d et al. I ll996t IColberg et al.lll999b lAltav et al.l 2006) . However, the orientations of dark matter halos 
are difficult to observe. One needs either X-ray observations for a sufficient number of groups and clusters 
of galaxies or reliable galaxy-group catalogs derived from the redshift surveys. The former are expensive 
and the latter are prone to a number of systematic errors. Despite these difficulties both a pproaches have 
been pursued. In most cases alignment signals out to, at least, 20 fe~ 1 Mpc have been found (lBinggelJ l982; 
lUlmer et alJl989l:l^s3ll989bHPlionisll994HChambers et al.l2000llHashimoto et al.l2007h . Instead of mea- 
suring the large-scale alignment based on groups and clusters of galaxies we suggest here to use directly 
the orientations of luminous galaxies and we quantify their alignment relative to large-scale structure. 

Galaxies are not oriented at random. Rather, the y have been found to show various forms of spatial 



alignment: between neighboring clusters of galaxies (Binggeli: 1982tlWestlll989aUPrionislll994l). between 
brighte st cluster galaxies (B CGs) and their parent clusters dCarter & Metcalfdll980l:1mnggelill982t[strabld 



1990; Hashim oto et"aT1l2 008). between the orientation of satellite galaxies and the orientation of the cluster 
dDekellll985t Plioni s et al. 1 120031) . and between the orientation of satellite galaxies and the orientation of the 



BCG (IStrubldl 19901) 



Observationally, these alignments are quantified either by the differential, P(9), or cumulative, P(9 < 
S mal ), probability distribution of the alignment angle 9 which is the angle between the major axis of a 
galaxy and the line connecting it to a neighboring galaxy. Also the mean values of those distributions, 
(9)(r p ), have been studied as a function of projected separation r p . With recent large redshift surveys, in 
particular the SDSS, it has become possible to determine the alignment using large and homogeneous sam- 
ples. Studies based on these surveys have focused mainly on the alignment of galaxies in groups. They 
revealed that satellite galaxies are preferentially distributed along the major axis of the central galaxies 
dBrainerdll2005l: lYang et alj|2006l lAzzaro et al]|2007t iFaltenbacher et alj|2007l), and satellite galaxies tend 
to be preferentially orien t ed toward the c entral galaxy ( Pereira & Kuhnll2005l : Agustsson & Brainerdf 2006: 



Faltenbac her et al.| [2007). Donos o et alj J200 6) analyzed a high-redshift sample (0.4 < z < 0.5) of lumi- 



nous red galaxies (LRGs) extracted from the SDSS, and found a clear signal of alignment between LRG 
major axes and the distribution of galaxies within 1.5/i _1 Mpc, indicating that the alignment effects observed 
in the local Universe were already present at z ~ 0.5. 

In this paper we propose two new statistics to quantify the alignment between galaxies and the large- 
scale structure. The first one we call the alignment two-point correlation function which is an extension of 
the traditional two-point correlation function. Basically, the correlation is measured as a function of pair 
separation and alignment angle. The second measure we call cos(29)-statistic and involves determining the 
mean cosine of twice the alignment angle for correlated pairs of given projected separation. We measure 
these statistics for SPSS galaxies as well as for sem i-analytic galaxies within the Millennium simulation 
(MS) dSpringeletal.ll2005H De Lucia & Blaizotl l2007l) . where the orientations of the semi-analytic galaxies 
are inferred from the orientation of their parent dark matter halos. The application of the new statistics 
reveals an alignment between red galaxies and large-scale structure out to 60 /i _1 Mpc. 

The paper is organized as follows. In § [2] we introduce the two statistics used here to quantify the 
alignment between galaxies and large-scale structure. The statistics are applied to the SDSS galaxy catalog 
in §[3] and compared to the results derived from the semi-analytic galaxy catalog of the MS in §[4] Finally, 
§|5] gives a short summary. 



2 METHODOLOGY 

For the subsequent analysis we introduce two new statistical measures to quantify the correlations between 
galaxy orientations and the cosmic density field. The first quantity we call the alignment (two-point) cor- 
relation function, w p (9p, r p ), where 9 V is the angle between the major axis of a galaxy and the connecting 
line to another one and r p is the projected separation between the two galaxies. This quant i ty is a two- 
dimensional extension of the traditional two-point projected correlation function. |Paz et al.1 d2008l) have 
used a related technique to analyze large-scale angular momentum alignments. The second quantity we 
refer to as the cos(29)-statistic . This statistic is closely related to similar quantities used in the context of 
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cosmic shear surveys. We compute both statistics only in projected space, i.e. orientations and pair separa- 
tions are two-dimensional vectors. We keep this restriction also in the second part of this analysis where we 
compare observational and numerical results. 

2.1 Alignment correlation function 

The two-point correlation function (2PCF) has long served as the primary way of quantifying the clustering 
properties of galaxies in redshift surveys (e.g. |Peebleslll980l) . It is defined as a function of pair separation 
by 

dP 12 =n 2 [l+Z(r)}dV 1 dV 2 , (1) 

where £(r) is the 2PCF, n the mean number density of galaxies, and dV\ and dV 2 are two infinitesimal 
volume elements centered at xi and X2 with separation r = X2 — xi. To be consistent with homogeneity 
and isotropy £ has been written as a function of the separation alone, that is, £(r). If £(r) ^ 0, then galaxies 
are said to be clustered. As the fundamental second-order statistic of the d ensity field, £ is sim ple to compute 
and provides a full statistical description for a Gaussian random field ( Bardeen et al.lll986l). It can al so be 
easily compared with the predictions of theoretical models (e.g. lAarseth et alll979tlDavis et al 1119851) . The 
amplitude of the correlation function on scales larger than a few Mpc provides a direct measure of the mass 
of the dark ma tter halos that host the galaxies through the halo mass-bias relation (e.g. lMo & Whitef l996: 
IJing et al.lll998l) . 

In galaxy redshift surveys, the 2PCF is measured in redshift space and expressed either as a function of 
redshift-space separation s, giving rise to a 2PCF of £(s), or as functions of separations perpendicular (r p ) 
and parallel (tt) to the line of sight, giving rise to £(r p , ir) with s 2 = tt 2 + r 2 . In many cases, the projected 
2PCF, w p (r p ), is the more useful quantity, as it does not suffer from redshift-space distortions, and is thus 
directly related to the real-space correlation function £(r). One can distinguish between two kinds of two- 
point correlation functions: two-point auto-correlation functions for which both members of a pair come 
from the same sample, and two-point cross-correlation functions (2PCCFs) for which the two members of 
a pair are from two different samples. In this paper we focus on the latter. In addition our analysis will be 
pursued in projected space, i.e. we focus on the projected 2PCCF w p {r p ) where r p denotes the projected 
separation of a galaxy pair. 

Given a sample of galaxies in question (Sample Q), a sample of reference galaxies (Sample R), and a 
random sample (Sample 1Z) that has the same selection function (i.e. distribution of redshifts and positions 
on the sky) as the reference sample, £(r p , tt) between Q and R can be estimated by 

N R QK(r p ,TT) 

where Nr and N-r are the number of galaxies contained in the reference and random samples, with 
Nh/Nr = 10 throughout this paper. QR(r p ,Tr) and QlZ(r p ,ir) are the counts of cross pairs between 
the indicated samples for a given separation perpendicular, r p , and parallel, tt, to the line-of-sight. With the 
measurement of £(r p , tt) in hand one can then estimate the projected cross-correlation function w p (r p ) by 
integrating £(r p , 7r) along the n direction: 

/+7T max 
£(r p ,7r)d7T = ^2^(r p ,TT t )ATTi, (3) 

where 7r max has to be sufficiently large to minimize the probability of erroneously excluding correlated 
pairs with line-of-sight separations larger than 7r max . 

Now we extend the definition of the 2PCCFs, so that they will be able to quantify the spatial alignment 
of galaxies. For each pair of galaxies with one member from Sample Q (the main galaxy) and the other 
from Sample R (the reference galaxy), we consider 9 p , the angle between the major axis of the main galaxy 
and the line connecting the two galaxies projected onto the sky. We include this angle as a second property 
of the pair, in addition to the pair separation. In this case, the correlation function is not only a function 
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of the projected separations, but also of 9 P . The estimator of Eq. is easily modified to account for the 
dependence on 9 p : 

g^r^.^m^-U (4) 

N R QK{9 p ,r p ,n) 

QR(9 p , r p , 7r) and Q1Z(9 P , r p , ir) are the counts of cross pairs between the indicated samples for given 9 p , 
r p and tt. The projected correlation function is found by integration along the line-of-sight. 

£(0 P , r p , n)dTT = ^2 €(Qp> r p> n i) AlT i ( 5 ) 

The traditional correlation function is just the average of the new correlation function over the full range 
of 9 p values. Taking symmetries into account the value of the angle ranges from zero (along the major 
axis of the main galaxy) to 90 degrees (perpendicular to the major axis). Thus, higher amplitudes of the 
new correlation functions at small values of 9 p indicate that the reference galaxies are more likely to be 
aligned along the major axis of the main galaxies. In contrast, higher amplitudes at larger angles indicate 
that the reference galaxies are more likely to be located along the minor axis of the main galaxies. This new 
statistic can be used for quantifying the alignment of galaxies and we refer to it as the alignment correlation 
function. 

On small scales (^ Mpc), this statistic can be used to confirm the alignment bet ween central and sa tellite 
galaxies in groups and cluster s, if Sample Q consists purely of central galaxies ("cf. lCarter & Metcalfe! 1980: 
Binggeli 1982: IStrublelll990h . More interestingly, the new statistic allows us to extend the alignment study 
to very large scales without worrying about selection effects, which are taken into account by comparison 
with the random sample (1Z). On large scales, this statistic can be used to quantify the alignment of the main 
galaxies (which may or may not be central galaxies) with respect to the large-scale structure of the Universe 
as probed by the large-scale distribution of the reference galaxies. 

As the final remark of this section we would like to mention that one can easily derive equivalent 
expressions for the alignment oMfo-correlation function - additionally it is straightforward to extend the 
formalism to three dimensional problems. 

2.2 The cos(26»)-statistic 

The cos(26>)-statistic gives the average of cos(2(?p) for all correlated pairs at a given projected separation. It 
will be referred as (cos(29 p ) cor ) (r p ) where the index cor emphasizes that the average is based on correlated 
pairs only. Again, 9 p indicates the angle between the major axis of a main galaxy and the line connecting it 
with a reference galaxy. More precisely, using the alignment correlation function, w p (9 p , r p ), we define 

, f00 , U s Io /2 cos{29 p )w p {9 pi r p )d9 p 

(cos(29 p ) cor )(r p ) = -y — — . (6) 

J w p (9 p ,r p )d9 p 

This statistic is constrained to values between —1 and 1, Values above and below indicate a preference 
for small (< 45°) and large (> 45°) angles, respectively. Values close to are expected for isotropy. An 
estimator for Eq. |6]is given by 



(cos(29 p ) cor )(r p ) = , , (7) 



V ( Q R e( r p> n ') \ Att- 
L>i \ Qn(r p ,TTi) J 



where QRe(r p , 7Tj) symbolizes the sum of cos(2# p ) for all cross pairs between main and reference samples 
(Q and R) within the given separation bins, r p and TTj. As before, QR(r p , 7Tj) and QlZ(r p , 7Tj) denote the 
number of cross pairs between the indicated samples for r p and TTj. 

Related statistics are used in weak lensing stud i es where they are ref erred to as ellipticity correlations 
(e.g. iMiralda-Escudd 1 1 99 ll: ICroft & Metzlerj bood iHeavens et~ai1l2000l) . In pa rticular we want to point 
out the similarity to the the intrinsic shear-density correlation function w g +(r p ) jMandelba um et al.l l2006: 
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iHirataet alJl2007h which measures the correlation between galaxy orientations and the large-scale den- 
sity distribution. L eaving the ellipticity weighting and 'responsivity' correction aside (cf. Eqs. 8 and 9 in 
iHirata et alJl2007h the following relation holds. 

w g + = w p {r p )(cos{20 p )) cor (r p ) (8) 

where w g + indicates the unweighted version of w g + . We do not weight by ellipticity since we here are 
solely interested in the spatial alignment between galaxies and the large-scale structure. 



3 ALIGNMENT OF SDSS GALAXIES 

In this section we apply the two new statistics to the publicly available data from the Sloan Digital Sky 
Survey DR6. Before doing that we describe some details of the survey and the galaxy sample construction. 
Also the determination of the galaxy orientations is reviewed. 

3.1 Galaxy sample construction 

The observational data used in this paper are taken from the SDSS which has been designed to obtain 
photometry of a quarter of the sky and spectra of nearly one million obj ects. Imaging is obtained in the u, 
g, r, i, z bands (Fukugita et al. 1996; S mith et alj|2002tllvezic et alj |2004). The details of the survey strategy 
can be found in lYork et al l d2000h and an overv iew of the data pipelines and products is provided in the 
Early Data Release paper fStoughton et alJl2002h. The SPSS has had seven ad ditional major data releases 
dAbazaiian et alj2003ll2004 120051: lAdelman-McCarthv et alj|2006ll2007ll2008h . 

The galaxy samples for this work are constructed from Sample dr6 of the New York University 
Value Added Galaxy Catalogue (NYU-VAGC), which is based on the SDSS DR6, publicly availa ble at 
http://sdss.physics.nyu.edu/vagc/ A detailed description thereof can be found in lBlanton et alJ d2005l) . Our 
sample consists of 430164 galaxies that are identified as galaxies from the Main sample (note that r-band 
magnitude has been corrected for foreground extinction), and are in the ranges of 0.01 < z < 0.4, —23 < 
Mo.i r < —17 and 14.5 < r < 17.6. Her e Mo.i r is the r-band absolute magnitude corrected to its z = 
0.1 value using the K — correction code of iBlanton et al.1 (l2003bl) and the luminosity evolution model of 



Blan ton et al.l ( 2003ah . We do not consider galaxies fainter than Mo.i r = — 17, because the volume covered 
by such faint samples is v ery small and t he results are subject to large errors as a result of cosmic variance 
(see for example Fig. 6 of lLi et ai1l2006l) . The faint apparent magnitude limit of 17.6 is chosen to yield an 
uniform galaxy sample that is complete over the entire area of the survey. This sample will serve as the 
reference sample (Sample R) for computing cross-correlation functions, as well as the parent sample for 
selecting different subsamples (Q). 

The large amount of data allows us to split the parent sample into various subsamples. Thus, we split all 
the galaxies into 4 subsamples according to their r-band absolute magnitudes, ranging from Mo.i r = —23 to 
Ma.i r = —17 with an interval of 1 magnitude. Further, we classi fy each galaxy with Ma.i r > —22 as either 



"red" or "blue" according to its g — r color. To this end, we follow lLi et al.l (120061) to fit the g — r distribution 
at fixed luminosity with a bi-Gaussian profile and use the mean of the two Gaussian centers as the color cut. 
In the highest luminosity bin, —23 < Mo.i r < —22, this color separation scheme is problematic, basically 
because the blue popul ation is extre mely sparse, thus we assign all galaxies within this magnitude bin to 
the red population (cf. lLi et ai]|2006l Fig. 9). These subsamples will be cross-correlated with the full parent 
sample, i.e. the reference sample. 

To compute the cross-correlation function between the main and reference samples, Q and R, one also 
needs to construct a random sample, 1Z, where galaxies with random coordinates are subject to the same 
selection effects as the reference sample. Since the reference sample is used as a tracer of the large-scale 
structure only the galaxy positions instead of their orientations are of interest. Accordingly, orientations are 
not considered when constructing the random samp le. A detailed account of the selection effects accom- 
panies the NYU-VA GC relea se (Blant on et al.ll2005l) which is the base for the construction of the random 
sample used here (cf. lLi et al]|2006l) . 



6 



A. Faltenbacher et al. 



3.2 Determination of galaxy orientations 

The orientation of galaxies is quantified by the position angle (PA) of the major axis of th eir r-band images, 
which is determined by the SDSS photometric pipeline called PHOTO dLupton et al.l200ll) . Photo provides 
three quantities for the PA of each galaxy: PAd e v, PA ea; p, and PA; so . The first two come from fitting two 
models to the two-dimensional image of the galaxy in each band: a pure de Vaucouleurs profile and a pure 
exponential profile, while the last one is given by measuring the shape parameters (centroid, major and 
minor axes, PA, and average radius) of th e 25 mag arcsec" 2 isop hote. Details on the photometric pipeline 
can be found in iLupton et al] d200ll) and IStoughton et ail d2002l) . In this paper we use the isophotal PA, 
PAi S0 . However, we found no significant variation in the alignment signals when adopting the alternative 
definitions for the PAs. 

To quantify the intrinsic scatter of the alignment signal we randomly shuffle the PAs among the galax- 
ies in the main sample and redo the alignment analysis. Since the shuffled orientations are not correlated 
with the large-scale structure one expects no systematic alignment signal. The la variance of 10 randomly 
shuffled samples can then be used to infer the intrinsic scatter. 

3.3 Alignment correlation function 

We now compute the alignment correlation function, w p (6 p , r p ), based on Eq.|5]with 7r max = 40 /i _1 Mpc. 
We probe projected separations up to 60 /i _1 Mpc within three angular bins. A less coarse angular binning 
is prohibited by Poisson noise at large separations. 

3.3.1 Alignment for red galaxies 

The first row of the upper panel of Fig. Q] shows the alignment correlation function for red galaxies in 
various magnitude bins. The red, green and blue solid lines display the results for the different angular bins, 
0° - 30°, 30° - 60° and 60° - 90°, respectively. The low angle bin, 0° - 30°, quantifies the correlation 
function along the orientation of the main galaxies. The 60° — 90° bin contains information about the 
reference galaxy abundance perpendicular to the orientation of the main galaxies. The long dashed black 
line represents the traditional correlation function and is a weighted average of the three bins. 

The second row of the upper panel of Fig. [T]shows the ratio of the alignment correlation function to the 
traditional one. The color code is the same as used in the row above. The shaded regions indicate the la 
variance between 10 random samples in which the orientations are shuffled at random among the galaxies. 
A signal well outside the shaded region is a significant detection of alignment or anti-alignment depending 
on whether the signal lies above or below one. 

The alignment correlation functions (colored lines) for red galaxies in the lowest luminosity bin (—20 < 
Mo.i r < —19) can hardly be distinguished from the traditional correlation function (dashed black line). 
However, the corresponding ratio, w p (6 p ,r p )/w p (r p ), shows a weak but systematic trend to anisotropy 
at scales larger than ~ 10 ft, _1 Mpc. At these separations the reference galaxies are preferentially located 
along the major axes of the main galaxies. Along the minor axes the reference galaxies are correspondingly 
under-abundant. This feature becomes more pronounced in higher luminosity bins. For galaxies brighter 
than Mo.i r = —20 a significant overabundance of reference galaxies is visible along the major axis. The 
signal reaches from 10 /i _1 kpc out to 60 ft," 1 Mpc which corresponds the entire range probed here. 

These plots indicate that the orientations of red galaxies are connected to the large-scale structure in 
which they are embedded. Any realistic galaxy formation model should be able to reproduce this interlacing 
between dimensions spanning ~ 3 orders of magnitude. 

3.3.2 Alignment for blue galaxies 

The lower panels of Fig. Q] display the corresponding results for blue galaxies. Otherwise the analysis is 
carried out in exactly the same manner as for the red galaxies. The panel for the highest magnitude bin, 
—23 < Mo.i r < —22, is left void because by definition there are no blue galaxies with these luminosities 
(see § 13.11 ). In the lower three magnitude bins we do not find any indication for alignment between the 
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-20.0<M. :r <-19.0 -21.0<M., r <-20.0 -22 0<M„, r <-21.0 -23.0<M„ r <-22.0 




r p [h *Mpe] r p [h 'Mpc] r p [h ! Mpc] r p [h ! Mpc] 




Fig. 1 Upper panels in each set: Solid lines display the projected SDSS alignment correla- 
tion function, w p (6 p ,r p ), between our reference sample and red (upper set) and blue (lower 
set) main galaxies in different intervals of r-band absolute magnitude. The colors correspond to 
three ranges in 9 p as indicated. The dashed black line shows the conditional correlation function 
averaged over angle. Lower panels in each set: Ratio between the angle dependent correlation 
functions and their mean. The color code is the same as used in the upper panels. The shaded 
regions plotted in red, green and blue colors indicate the 1 — a variance based on 10 samples in 
which the position angles are shuffled at random among the galaxies. By definition (§ 13. Il l there 
are no blue galaxies in the brightest luminosity bin so these panels are kept void. 



galaxies and the reference galaxy distribution. The differences seen between red and blue galaxies suggest 
that the orientations of red and blue galaxies are determined by distinct physical processes. 

3.4 cos(26»)-statistic 

In this section we compute the cos(26>) statistic based on the same subdivision of the parent SDSS sample 
by color and r-band absolute magnitude as used above for the alignment correlation function. In a second 
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Fig. 2 cos(26> p ) statistic for SDSS data. The results for red and blue main galaxies are displayed 
by red and blue lines. The galaxies are subdivided into four r-band absolute magnitude bins as 
indicated. The shaded regions in red and blue show the ler variance between 10 random samples 
in which the position angles are shuffled at random among the main galaxies. The reference 
sample is the same for all these measurements and comprises all galaxies with r-band absolute 
magnitudes in —23 < Mi.i r < —17. No results are shown for blue galaxies in the highest 
luminosity since, by definition, they are red galaxies. 



step we focus on the alignment signal for group central galaxies. This is to facilitate a comparison with 
simulation results discussed in the second part of this paper. 

3.4.1 Color and luminosity dependence 

Fig.|2]displays the cos(2# p )-statistic as discussed in § 12. 21 The results for red and blue galaxies are displayed 
by red and blue lines. The shaded regions show the la variance between the 10 random samples. The signal 
for the red galaxies shows a strong dependence on luminosity. There is basically no signal for red galaxies 
in the lowest luminosity bin. For r-band absolute magnitudes, ranging from Mo.i r = —20 to Mo.i r = —23 
the signal systematically lies above the shaded region, indicating that it can not be explained by the intrinsic 
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Fig. 3 cos(2#p) statistic for SDSS central galaxies. See the text for our definition of 'central'. 
Results for red and blue galaxies are displayed by red and blue lines. All symbols are the same 
as used in Fig. [2] 



scatter inherent to the data. The signal is significant on all scales probed (< 60 ft. _1 Mpc) and becomes more 
pronounced with increasing luminosity. Blue galaxies do not show any indication for alignment with the 
reference galaxy distribution. Fig.[2]suggests that red galaxies with L > tend to be aligned with large- 
scale structure out to at least 60 ft, _1 Mpc. This result is very similar to that from the alignment correlation 
function discussed in § 13.31 It supports the picture of a physically distinct origin for the orientations of red 
and blue galaxies. 



3.4.2 Alignment of group central galaxies 

Alignment between clusters of galaxies has been detected out to large scales ~ 100 h~ 1 Mpc) (e.g jBin ggeli 
1 9821 IPlionisll 1 994 lHashimoto et al.ll2007h . In addition there exists strong evidence that the orientations of 
central galaxies are aligned with the orientations of their parent groups or clusters. Now the question arises 
whether an exclusion of non-central galaxies may enhance the alignment signal measured before. To this 
end we redo the above analysis excluding non-central galaxies. This approach also facilitates a subsequent 
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comparison with N-body results, since there only the orientations of central galaxies are available (see 

Central galaxies are found in the following way. For each main galaxy, we use its r-band absolute 
magnitude to ca lculate a ha l o mass , according to the relation between central galaxy lu minosity and hal o 
mass as given bv lYang et al](l2005al) . We then estimate a 'virial' radius using the model of lEke et al.l d200ll) . 
We find all the companions around the galaxy within this radius and a line of sight velocity separation 
of ±3000 kms" 1 and compare the luminosity of the galaxy to that of the companions. If the galaxy is 
brighter than any of its companions and it is not inside the virial radius of a farther, brighter galaxy it is 
classified as a 'central' galaxy. Otherwise, it is a 'non-central' galaxy. This method will inevitably mis- 
judge a certain fractio n of galaxies to be non-central due to the appearance of more luminous interlopers 
(cf. lYangetal.ll2005bh . However, it should be sufficient to estimate the signal for central galaxies. This can 
then be compared to the signal for the complete sample as well as to simulations. 

Fig. [3] displays the alignment signals based on central galaxies only. As before we do not find any 
alignment signal for blue galaxies. The upper two luminosity bins show only marginal differences when 
compared to the complete sample. This can be explained by the fact that most bright galaxies are central. 
For the lower two luminosity bins the alignment signal increases by a few per cent which indicates a certain 
reduction of the alignment signal for the complete sample which includes non-central galaxies (Fig. [2j». We 
come to the conclusion that the alignment signal is not strongly enhanced by the exclusion of non-central 
galaxies which is due in part to the fact that luminous galaxies usually are central. 

4 COMPARISON TO THE ALIGNMENT OF MS GALAXIES 

In the second part of this study we apply our new analysis tools (alignment correlation function and cos(2#)- 
statistic) to semi-analytic galaxies in the MS. Orientations have to be assigned to the model galaxies and 
we investigate different approaches, comparing the resulting alignment signals with the observed ones. 
Since we did not detect any alignment signal for blue galaxies we restrict our orientation assignment to red 
model galaxies. This restriction may also be justified by a slightly different kind of reasoning. Namely, red 
luminous galaxies are predominantly elliptical galaxies. Those are thought to form via merging of proto- 
galaxies which is basically a collisionless process. Thus N-body simulations may be able to recover some 
characteristic features of red galaxies, in particular their orientation. Blue galaxies, on the other side, are 
mainly spiral galaxies with their properties heavily dependent on baryonic, i.e. collisional, physics. Thus 
their properties may be only poorly recovered by N-body simulations. This picture may change at high 
redshifts. At that time, also the orientation of many blue galaxies could be a product of gas-rich mergers. 

In the following two paragraphs we list some details of the MS and describe how we estimate orienta- 
tions for the model galaxies. Subsequently we redo the analysis already carried out on SDSS galaxy samples 
and we compare the results. 

4.1 The simulation and the galaxy sample 

The Millennium Simulation (Springel 2005) adopted concordance values for the parameters of a flat A cold 
dark matter (ACDM) cosmological model, 57d m = 0.205 and fib = 0.045 for the current densities in CDM 
and baryons, h = 0.73 for the present dimensionless value of the Hubble constant, erg ~ 0.9 for the rms 
linear mass fluctuation in a sphere of radius 8 h~ 1 Mpc Mpc extrapolated to z = 0, and n = 1 for the 
slope of the primordial fluctuation spectrum. The simulation followed 2160 3 dark matter particles from z= 
127 to the present day within a cubic region 500 /i _1 Mpc on a side resulting in individual particle masses 
of 8.6 x 10 s h~ 1 MQ. The gravitational force had a Plummer -equivalent comoving softening of 5 /i _1 kpc. 
The Tree-PM iV-body code GADGET2 (Springel etaL||2005j) was used to carry out the simulation and the 
full data were stored 64 times spaced approximately equally in the logarithm of the expansion factor. This 
information makes it possible to construct trees that store detailed assembly histories for each dark matter 
halo present at z= 0. 

The halos are found by a two-step procedure. In the first step all collapsed halos with at least 20 
particles are identified using a friends-of-friends (FoF) group-finder with linking parameter b = 0.2. These 
objects will be referred to as FoF-halos. Then post-processing with the substructure algorithm SUBFIND 



Alignment between galaxies and large-scale structure 



11 



M r 


[-20,-19] 


[-21, -20] 


[-22, -21] 


[-23, -22] 


Ahalo 


16027 


98618 


167513 


53065 


<AU) 


489 


1672 


4707 


18662 


(iVcen) 


31 


83 


164 


278 



Table 1 Particle number statistics for host halos 



(Springeletafl[2001) subdivides each FoF-halo into a set of self-bound sub-halos. Here we only consider 
the main sub-halo which is the sub-halo with the most massive progenitor among all sub-halos belonging 
to the same FoF-halo. We refer to this sub-halo as the parent sub-halo associated with the central galaxy. 

Based on the assembly histories individual halos are populated with semi-analytic galaxies for which 
many 'observable' quantities are generated. For a detailed des cript ion of the construct i on of the semi- 
analytic galaxy catalog we refer the reader to lCroton et alj d2006l) and lDe Lucia & Blaizoj J2007). Here we 
use the DeLucia2006a_SDSS2MASS catalog flhttp : / /www . g-vo . org/MyMillennium2 /[ ) which 
provides synthetic magnitudes through SDSS filters, alternatively to classify the model galaxies as well as 
to tag their parent sub-halos. The use of semi-analytic galaxies facilitates comparison to the observations 
and is also an elegant way to characterize the assembly history of the host halo by a few numbers. 

Our standard reference sample comprises a random subset (1103640 galaxies) of all semi-analytic 
galaxies with 7--band absolute magnitudes of —23 < M r < —17 (in total 1 1027979 galaxies) at z = 0. This 
standard reference sample is used for all the subsequent analysis based on MS galaxies. The parent galaxy 
sample comprises all red, central galaxies within the same magnitude range as the reference galaxies ( 
—23 < M r < —17). A galaxy is red if g — r > 0.7, where g and r are synthetic magnitudes in the 
corresponding SDSS bands. A galaxy is central if it is the dominant galaxy in a given FoF-halo. Expressed 
in the MS terminology it is the (only) type galaxy within this FoF-halo. This definition comes close to 
the selection we have employed to identify central galaxies in the observational analysis above. In analogy 
to the SDSS sample the MS galaxies are split into four r-band magnitude bins with equal intervals of 
one magnitude. Each model galaxy is associated with a dark matter halo which will be used to determine 
the orientation of the galaxy. For the four magnitude bins Table Q] lists: the number of halos (A^aio); the 
average number of particles belonging to these halos, (iV su b } ; and the average number of particles belonging 
to the central part of the halos, (N ccn ). The determination of the number of central particles is related to 
the computation of the central orientations and will be discussed below. We conclude this paragraph by 
emphasizing that orientations are only determined for red, central galaxies, because only for those galaxies 
do we expect a tight correlation between the orientations of galaxy and halo. 



4.2 The orientations of semi-analytic galaxies 

As a proxy for the major axis of the ce ntral galaxy we adopt the major axis of the projected moment of 
inertia tensor of the parent sub-halo (cf. Agustsson & Brainerd 200(1 iKang et aT]|2007l Faltenb acher et alj 



2008; iKnebe et al.ll2008t lOkumura et al.1 12008). Since real galaxies are seen in projection, we project the 
dark matter distribution before diagonalizing the moment of inertia tensor. We determine the orientations 
of each central galaxy in two alternative ways. We use the dark matter distribution of the parent sub-halo, 
(i.e. the main sub-halo of the FoF-halo) or we only use its central part to approximate the orientation of the 
central galaxy. The former is d one to compare the results to earlier work w hile the latter better mimics the 
true physical circumstance (cf . iFaltenbacher et al.ll200 8; Knebe et al. 2008). We will refer to the former as 
halo orientation and the latter as central or galaxy orientation. The central orientation is determined in an 
iterative way. We begin by computing a number of central particles iV cen using the following formula 

v _,- x^6 N sub > 4900 
" xN™ for N sub < 4900 W 

where N su b is the number of particles within the parent sub-halo. This formula is a two-power-law approx- 
imation to the observed relation between the luminosity of a central galaxy and its parent halo mass (e.g. 
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Fig. 4 Probability distribution of the misalignment angle between halo and central (galaxy) 
orientations for a series of magnitude bins as indicated. Errors represent Poisson uncertainties. 
The small vertical lines on the bottom axis denote th e mean values of the distributions. The 
dashed line display the result of lOkumura et alj J2008) and its mean we show as black dashed 
mark on the bottom axis. 



ICoorav & Milosavrjevicll2005l) . It reflects the assumption that the distribution of a galaxy's stellar com- 
ponent can be represented by the central matter distribution of its simulated dark matter halo. The mass 
representing the galaxy then has to follow the observed scaling laws. Now, let be the radius which 
encloses the N cen . In a first iteration step the moment of inertia is computed within the r^t. The resulting 
axis ratios are used to cut out a central ellipsoidal where the length of the intermediate axis is fixed to ri n a . 
Based on this new subset of particles the moment of inertia is computed anew. The iteration is continued 
until the central orientation converges. Fixing the intermediate axis results in approximately unchanging 
particle numbers within consecutive ellipsoids. Up to this point the computations were done in three dimen- 
sions. In a last step all particles within the final ellipsoid are projected onto the plane of the sky. The inertia 
tensor of this 2-D distribution determines the central orientation. 

In Fig.|4]the difference between halo and central (galaxy) orientations are illustrated by the probability 
distribution function (PDF) of the misalignment angle between them. The mean values of the misalignment 
for each magnitude bin are indicated by the little vertical marks on the bottom axis. As indicated in TableQ] 
on average there are only 3 1 particles used to determine the central orientations within the lowest luminosity 
bin. Thus, Poisson noise adds an additional misalignment component which contributes to the flattening 
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Fig.5 Upper panels: Projected MS alignment correlation function, w p (9 p , r p ), for red central 
galaxies in different intervals of r-band absolute magnitude using the standard reference sam- 
ple. The colors correspond to the sectors as indicated. The solid and dashed lines are based 
on central and halo orientations. The dotted line represents the traditional projected two-point 
cross-correlation function, w p (r p ). Lower panels: Ratio between the alignment and traditional 
correlation functions. The color code is the same as used above. Red and blue shaded regions 
indicate the 1 — a variance based on 10 samples in which the central orientations are shuffled at 
random among the galaxies. 



of the PDF at small angles. In fact calculating the misalignment using central volumes twice as large as 
given by Eq. [9] substantially reduces the flattening in the lowest luminosity bin. However, since we aim to 
probe the volume actually occupied by the stellar component of the galaxies we adhere to Eq. [9] keeping 
in mind that the central orientations for the low luminosity galaxies are poorly resolved. Such effects are 
less important for the higher luminosity bins. The two highest luminosity bins use 160 an d 270 particles on 
average which should result in a robust determi nation of the ori e ntatio ns (cf. |jinal2002l) . For comparison 
we also display the Gaussian PDF suggested bv lOkumura et alj J2008), who found a mean misalignment 
angle of ~ 25° between the halo and the galaxy orientation. Comparing the orientations of cent ral galaxies 
and th eir host halos a similar misalignment ha s been found bvlAgustsson & Brainerdl d2006l) and|Kang et al. 
d2007h . Somewhat lower values are quoted in lWang et alj d2008l) . The misalignment weakens slightly with 
increasing luminosity of the central galaxy, though this may be due to the resolution effects discussed above. 

4.3 Alignment correlation function 

With these orientations and the model magnitudes the alignment, correlation function can be computed for 
the galaxies in the MS. We will adopt the same luminosity bins used for our analysis of the SDSS data, 
simplifying the comparison between observations and simulations. For the computation of the alignment 
angle we apply the distant observer approximation. To that end the z-direction in the simulation is chosen 
to be parallel to the line of sight. To mimic the observational approach the maximum projected separation 
of pairs along the line of sight is 7r max = 40 ft, _1 Mpc. 

Fig.|5]shows the alignment correlation function for four luminosity bins which can be compared to the 
upper panel in Fig. Q] Also the luminosity range of the reference sample is chosen to be the same as in our 
observational analysis. At this point we would like to emphasize that the MS main sample only includes 
central galaxies. This causes the pronounced bend at a few hundred /i _1 kpc corresponding roughly to virial 
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radii of the host halos. The general behavior, however, is very similar to that seen in Fig.Q] The correlation 
function is enhanced along the major axes of the red galaxies (0° < 8 p < 30°) and reduced perpendicular 
to them. For galaxies brighter than M r < — 20 a systematic alignment/anti-alignment signal is seen over 
the entire separation range probed (out to 60 /i _1 Mpc). And the strength of the alignment effect increases 
with the luminosity of the galaxies. 

In addition we find that, except for small (intra-halo) scales, the alignment signal based on the halo 
orientations is more pronounced than that based on the central (galaxy) orientations. This can be explained 
by the misalignment between halo and central orientations discussed in § 14.21 A comparison with the upper 
panel of Fig.Q] re veals that the alignment based on halo orientations overestimates the observational signal 
whereas that based on the central orientations reproduces the amplitudes obtained for SDSS galaxies quite 
well. 

So far we have restricted the computation of the alignment correlation to 3 angular bins, a relatively 
coarse segmentation. This facilitates the direct comparison with the SDSS results, derived above, where a 
less coarse binning would result in excessive Poisson noise. However the MS provides much better statistics, 
making possible the computation of a 2D alignment correlation map as displayed in Fig. [6] This shows the 
alignment cross-correlation between a sample consisting of red (g — r > 0.7), central galaxies with r-band 
luminosities between -21 and -22 and our standard reference sample. The left and right panels are based 
on halo and central orientations, respectively. The panels in a column display consecutive blowups of the 
same map. The correlation amplitude can be read off using the adjacent color bars. Note, the range changes 
from top to bottom. The ratios of minor to major axes for log w p (x,y) = 1.5, 1 and 0.5 are 0.57, 0.61 
and 0.42 for the halo orientations and 0.75, 0.80, 0.75 for the central orientations. In agreement with the 
results above, the alignment correlation based on halo orientations is more anisotropic. It is striking that the 
flattening of the counts is almost constant over the range of scales we can measure. 

The main conclusion to be drawn from Figs. [5] and [6] is that the halos of red luminous galaxies are 
aligned with the surrounding galaxy distribution out to at least 60 h~ 1 Mpc. If the orientations of the galax- 
ies are approximated by the orientation of the halos as a whole, the alignment signal is strong. The mis- 
alignment between the halo and central region, as discussed in § 14.21 causes a weakening of the correlation. 
The detection of the alignment between red galaxies and the large-scale structure provides a tool to test 
predictions for the alignment of galaxies with their halos. 

4.4 cos(26»)-statistic 

To facilitate the comparison of observed and modeled cos(26>)-statistics we follow the observational ap- 
proach discussed in § 13.31 i.e we use the equivalent luminosity and color cuts. Again the analysis is carried 
out in projected space where the z-axis is chosen to lie parallel to the line-of-sight and the distant observer 
approximation is applied. The maximum line-of-sight separation between correlated pairs is restricted to 
40 /i _1 Mpc. Only the signal based on red (g — r > 0.7), central galaxies is compared to the corresponding 
SDSS results. The standard MS reference sample is used. 

In Fig.|7]the results for central, red SDSS galaxies are copied from Fig. [3] including the shaded region 
reflecting the intrinsic scatter in the observational data. The black lines display the corresponding signals 
derived from the MS. Thick lines represent the results based on central orientations, whereas thin lines 
are based on halo orientations. In the lowest luminosity bin the halo orientations result in a substantial 
alignment signal whereas the signal based on the central orientations and the SDSS data are in agreement 
with no detected alignment signal. For the next two luminosity bins the signal based on central orientations 
agrees well with the SDSS data. The mismatch for the largest separation bin may be associated with zero- 
crossing of the correlation function, because the response of the cos(26>)-statistic is very sensitive near this 
point. In the highest luminosity bin there appears a disagreement between the MS signal based on central 
orientations for separations < 5 h~ 1 Mpc and the SDSS measurements. The agreement is better with the 
MS results based on the halo as a whole. On larger scales 5 /i _1 Mpc < r p < 60 /i _1 Mpc again there is 
good agreement with the signal based on central orientations whereas that based on overall halo orientations 
is an overestimate by 50-100 per cent. 

We are led to the following conclusions: 1) the alignment signal depends on the magnitude of the central 
galaxy and/or the mass of the host halo. 2) The alignment signal based on the orientation of the halo exceeds 
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Fig. 6 2D alignment correlation function based on the DeLucia2006a_SDSS2MASS semi- 
analytic galaxy catalog for the MS. The main sample consists of red, central galaxies with abso- 
lute r-band magnitudes between —21 and —22. The orientations of the galaxies are aligned with 
the x-axis. Left and right panels display results based on halo and central orientations. Panels 
from top to bottom show the same distribution for different scales. The color code indicates the 
correlation amplitude, its range changes from top to bottom. Pixels with values below and above 
the indicated range are displayed black and yellow, respectively. The ratios of minor to major 
axes at log w p (x, y) = 1.5, 1 and 0.5 are 0.57, 0.61 and 0.42 for the halo and 0.75, 0.80, 0.75 for 
the central orientations. 



the observed signal, in most cases, whereas the central orientations generally result in good agreement 
between observations and simulations. 3) The alignment for red central galaxies with Mo.i r < —20 persists 
to the largest scales probed (~ 60 ft, Mpc). 
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Fig. 7 Comparison of the cos(2#) signal for red SDSS galaxies and semi-analytic galaxies. The 
red lines and symbols are copied from Fig. [3] The black lines display the corresponding results 
derived from semi-analytic galaxy catalogs. Thin lines are based on the halo-orientations and 
thick lines represent results based model central-orientations. 



5 SUMMARY 



In the first part of this study we use the SDSS DR6 galaxy catalog to quantify the alignment between 
galaxies and the large-scale structure. We develop two new statistics, namely the alignment correlation 
function and the cos(29)-statistic. The former is a two-dimensional extension of the traditional two-point 
correlation function taking into account the orientations of the main galaxies. The latter is related to the 
ellipticity correlation function used to analyze cosmic shear observations. In both cases we represent the 
large-scale structure by a reference galaxy sample with absolute r-band magnitudes between —17 and —23. 

The alignment correlation function is defined on the 2D plane of the sky. In our SDSS analysis we use 
a polar coordinate system and probe the correlation function along the galaxy major axes focusing on the 
0° < 6 P < 30° sector, where 6 p is the angle between the major axis and the connecting line to a reference 
galaxy. The abundance of reference galaxies perpendicular to the major axes is represented by the counts 
within the 60° < 9 p < 90° sector. The main results revealed by the alignment correlation function are: 1) 
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The major axes of red galaxies with Mo.i r < —20 are significantly aligned with large-scale structure; 2) 
The signal is visible out to 60 /i _1 Mpc (the range probed here); 3) Independent of luminosity blue galaxies 
are not aligned with large-scale structure. 

The cos(26>)-statistic also probes the anisotropy of the reference galaxies around main galax i es and 
is r elated to the galaxy — intrinsic shear correlation function as discussed in iMandelbaum et all (12006b 
and Iffirata et al.l ( 12007 ). In contrast to the alignment correlation function it does not require angular bin- 



ning. It also reveals significant alignment between red luminous galaxies and large-scale structure out to 
60 h~ 1 Mpc but no alignment is detected for blue galaxies. The restriction to central galaxies results in a 
slight enhancement of the alignment signal for red galaxies, in particular, at lower luminosities and at large 
separations. However, the differences are small, indicating littl e reduction of the sig nal by satellite galaxies 
which may have orientations altered by the local tidal field (cf. Per eira & Kuhnll2005h . 

The second part of this study is based on a Millennium Simulation semi-analytic galaxy catalog where 
we have assigned an orientation to each red, luminous and central galaxy based on the orientation of the 
mass distribution of the inner halo. As an alternative identification of the orientation of the central galaxy 
we also used the orientation of the host halo itself. 

The mean projected misalignment angle between a halo as a whole and its central (galaxy) region is 
~ 25°. This misalignment decreases slightly with increasing luminosity of the central galaxy and/or its 
host halo mass. This behavior is particularly notable, since the relative size of the central region (used to 
determine the central orientation) decreases with halo mass (cf., Eq. |9}. 

Based on the orientations and the luminosities of the semi-analytic red galaxies we compute the align- 
ment correlation function in the same manner as for the SDSS galaxy catalog. We find that the alignment 
signal based on the halo orientations overestimates that found for SDSS. With the central orientation, how- 
ever, quite good agreement between the modeled galaxy sample and the observations has been achieved. 
In other words, a misalignment between halo and central orientations is crucial to matching the observed 
amplitude of the alignment signal. Our approach to compute the major axis of the central elliptical galaxy 
based on the central dark matter distribution is physically motivated and results in the right amount of mis- 
alignment. However, any other process which leads to equal misalignment will be in agreement with our 
measurements as well. 

The large data volume of the MS allows us to generate a two-dimensional map of the alignment cor- 
relation function w p (x, y), i.e. we can adopt a very fine two dimensional binning and still get a reasonable 
number of pair-counts for each bin. The maps based on halo orientation show a substantially flattened 
galaxy distribution around red main galaxies out to at least 60 7i -1 Mpc. Using central orientations the flat- 
tening is reduced. In particular we have shown a map based on red main galaxies with ? -band magnitudes 
between —22 and —21 using the standard reference sample. There the ratios of the minor to major axes at 
\ogw p (x, y) — 1.5, 1 and 0.5 are 0.57, 0.61 and 0.42 for halo orientations and 0.75, 0.80, 0.75 for central 
orientations. It is interesting that the flattening is almost independent of scale all the way from 1 h~ 1 Mpc 
to 60 h^Mpc. 
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